Two time scales and FDT violation in a Finite Dimensional Model for Structural 

Glasses 
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We study the breakdown of fluctuation-dissipation relations between time dependent density- 
density correlations and associated responses following a quench in chemical potential in the Frus- 
trated Ising Lattice Gas. The corresponding slow dynamics is characterized by two well separated 
time scales which are characterized by a constant value of the fluctuation-dissipation ratio. This 
result is particularly relevant taking into account that activated processes dominate the long time 
dynamics of the system. 



In recent years considerable progress has been achieved 
in the theoretical description of the glassy state of mat- 
ter. A scenario for the observed slow dynamics of glass 
forming materials has emerged through detailed analysis 
of mean field (MF) spin glass models pj. The equa- 
tions describing the off-equilibrium dynamics of these MF 
spin glasses simplify, above the transition, to the single 
equation for the Mode Coupling Theory for supercooled 
liquids 0|. These approaches have been successful in 
explaining history dependence or aging effects and the 
nature of the two characteristic relaxations in glasses, 
the short time or /3-relaxation and the structural long 
time a-relaxation. In the a-relaxation the system falls 
out of local equilibrium, as a consequence Fluctuation- 
Dissipation Theorem (FDT) breaks down and can be re- 
placed by the more general relation 



R(t, t w ) 



X(t,t w ) dC{t,t w ) 



T 



dt. 



(i) 



where C(t,t w ) is a two times correlation function and 
R(t,t w ) with t>t w is the associated response. T is the 
heat bath temperature and X{t,t w ) is a function that 
measures the departure from FDT: at equilibrium X = 1 
and the usual FDT is recovered, while in the out of equi- 
librium regime X < 1. In MF approximation || the 
function X, called "fluctuation-dissipation ratio" (FDR), 
turns out to depend on both times only through C(t, t w ). 
Moreover in MF models of glasses A is a constant (when 
different from 1). This scenario reflects the existence of 
only two well separated time scales, the equilibrium or 
FDT scale and a longer one where the system is out of 
equilibrium. The FDR has been interpreted as an effec- 
tive temperature and it has been demonstrated that it is 
exactly the temperature that a thermometer would mea- 
sure if it would be coupled to the slowly relaxing modes 
of the system M . Recently the first experimental deter- 
mination of the FDR has been done in glycerol pi. 



This scenario, while appealing, is essentially based on 
an analogy between the physics of some MF spin glasses 
and the behavior of real finite structural glasses (being 
the formal connection valid only in the high temperature 
region). At this point it seems crucial to test the link 
between MF theories and realistic models in the glassy 
phase. In particular, we still do not know which will be 
the role of activated processes in realistic models. Acti- 
vated processes are absent in completely connected mod- 
els in the thermodynamic limit while on the contrary 
they dominate the relaxation dynamics below the glass 
transition temperature T g in real glasses. 

The goal of the present letter is to go beyond the MF 
like description of structural glasses by considering a fi- 
nite dimensional model, the frustrated Ising lattice gas 
(FILG) d,EJ, wn i cn presents most of the relevant fea- 
tures of glass forming materials, in particular activated 
processes at low temperatures. Here we analyze the vio- 
lation of FDT in the FILG in three dimensions through 
Monte Carlo simulations. The (very precise) results con- 
firm the qualitative scenario of MF models of a constant 
FDR with large separation of time scales and set the 
stage for a detailed investigation of activated processes 
in realistic models of glasses. 

The FILG is defined by the Hamiltonian: 



H = - J } (e i: jO-iCTj - l)ninj - 



(2) 



At each site of the lattice there are two different dy- 
namical variables: local density (occupation) variables 
rii = 0, 1 (i = 1 . . . N) and internal degrees of freedom, 
(Tj = ±1. The usually complex spatial structure of the 
molecules of glass forming liquids, which can assume sev- 
eral spatial orientations, is in part responsible for the 
geometric constraints on their mobility. Here we are in 
the simplest case of two possible orientations, and the 
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steric effects imposed on a particle by its neighbors are 
felt as restrictions on its orientation due to the quenched 



random variables Si 



±1. The first term of the Hamil- 



tonian implies that when J — > oo any frustrated loop in 
the lattice will have at least one hole and then the density 
will be p < 1, preventing the system from reaching the 
close packed configuration. The system will then present 
"geometric frustration" . Finally, p represents a chemical 
potential ruling the system density (at fixed volume). 

The system presents a slow "aging" dynamics after 
a quench from a small value of p characteristic of the 
liquid phase to a large p corresponding to the glassy 
phase 0, what is equivalent to a sudden compression. 
In the present numerical experiments we always let the 
system evolve after a quench in p with J and T fixed. 
The origin of the times is set on the quench time. Af- 
ter the quench the density slowly relaxes up to a critical 
value near p c * 0.675 (further details will be given in ||). 
After a waiting time t w we fix the density to the value 
P = p(tw) (for technical reasons explained below) and a 
small random perturbation (/ij = ±1) is applied M: 



H'(t) = H(t) - e{t)J2»iMt) 



(3) 



In all our numerical experiments the field is switched 
on at time t w and kept fixed for later times, that is 
e(t) = e8(t — t w ). Then we measure the density-density 
autocorrelation function 



C(t,t w ) = — ^2 (ni{t)ni(t w )) 



(4) 



where ( • } and ~ are the averages over thermal histories 
and disorder realizations |IFJ] . At the same time we mea- 
sure the associated response function integrated over the 
time and divided by the perturbing field intensity, which 
defines the off-equilibrium compressibility 

K (t,t w ) = - [ R(t, s)e(s)ds = [ R(t, s)ds , (5) 



(6) 



where, as usual, the response is defined as 



R(t, t') 



Np 



E 



d(nj(t)) 
de(t') 



In the large times limit (t, t w — > oo), X(t, t w ) depends on 
both times only through the correlation C(t,t w ). Then 
integrating Eq. ([!]) from t w to t we obtain a useful relation 
linking the correlation and the compressibility in the out 
of equilibrium regime 



(7) 



Tn(t,t w ) = / X{C)dC 
JC(t,t w ) 



In our case the perturbing term in the Hamiltonian, 
shown in Eq. (0) , gives to the integrated response the fol- 
lowing form 



(8) 



where [ • ] av is the average over the random pi realiza- 
tions []l0|. The second term can be ignored because the 
Pi are random and completely uncorrelated from the con- 
figuration at time t w . 

Performing a parametric plot of the compressibility (or 
the integrated response) versus the correlation is a useful 
way of getting information about the different dynamical 
regimes present in the model and in particular the time 
scales structure of the system. In fact from Eq. (Q) it is 
easy to see that, plotting Tre(i, t w ) vs. C(t, t w ), the FDR 
can be simply obtained as minus the derivative of the 
curve, i.e. 



X{C) 



d[TK(t,t v 



(9) 



C(t,t w )=C 



This is the key relation used in order to extract the FDR. 



There is already a considerable literature on this kind 
of analysis in systems with and without quench disor- 
der [jll]]. Particularly relevant to the present discussion 
are references |l2[, where a constant FDR was found 
in model glasses with interactions of the Lennard- Jones 
type and in a purely kinetic lattice gas. The FILG has 
the advantage of being a Hamiltonian lattice model with 
short range interactions and, in this sense, it is more real- 
istic than purely kinetic models and more accessible an- 
alytically and computationally than Lennard- Jones sys- 
tems. Moreover, it is a valid on-lattice model for struc- 
tural glasses and it may be simple enough to apply sta- 
tistical mechanics techniques. 

We have simulated the FILG in 3D for linear sizes 
L = 30 and 60. Fixing the coupling constant J = 1 
and a temperature T = 0.1 the system presents a glass 
transition around p = 0.5 ||. 

In all our numerical experiments we have prepared the 
system in an initial state with low density, characteristic 
of the liquid phase and then, at time zero, we have per- 
formed a sudden quench in p to a value deep in the glass 
phase (the data presented here refer to p = 1). As al- 
ready explained above, after a time t w a perturbation in 
the form of a random small chemical potential has been 
applied and the density-density autocorrelation, Eq.(^), 
and the corresponding integrated response, Eq.(|^), have 
been recorded. At time t w the density has been fixed to 
the value p = p(t w ) for the following reason: the perturb- 
ing term in the Hamiltonian [see Eq.(|^)] favors roughly 
half of the sites (those with pi = 1) and the remaining 
half to be empty. When the perturbation is switched 
on, the first half starts to be more filled than the sec- 
ond one and this increases the integrated response, as it 
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should. On very long times, however, because the den- 
sity reached at time t w is still a bit below the asymptotic 
value, the number of particles continues growing and new 
particles are added with higher probability on the sites 
with [ii = — 1, which are more empty. Then the response 
becomes negative because of systematic errors. We avoid 
the negative responses fixing the density at time t w . In 
the limit t w — > oo we recover the right behavior in any 
case, however with our choice the extrapolation is safer. 

In all our simulations we have verified that the strength 
of the perturbation e = 0.02 is small enough in order to 
be in the linear response regime. We have also checked 
that the system thermalizes at temperatures T > 0.3 (be- 
ing always fi = 1) and satisfies FDT. Further details will 
be given in ||. 




0.7 0.75 0.8 0.85 0.9 0.95 1 
C(t,tJ 

FIG. 1. The plot of the integrated response times the tem- 
perature versus the correlation gives clear evidence for the 
existence of two well separated time scales, together with a 
constant FDR in both regimes. The dashed line is Tk = 1— C. 
The errors are of the order of the symbols and have been es- 
timated from sample to sample fluctuations. 

In Fig. |]we show the main result of this letter: a para- 
metric plot of the integrated response versus the density 
autocorrelation for different waiting times. The behavior 
of the curves is exactly the one predicted by MF the- 
ories for a glass former. Two distinct regimes can be 
perfectly recognized. In the FDT or quasi-equilibrium 
regime, t — t w <C t w (/3-relaxation), the points lie on the 
straight line given by 

T K {t,t w ) =1-C{t,t w ) . (10) 

This first time regime corresponds then to a FDR 
X(t,t w ) = 1 independent of t and t w . In this regime the 
system is in quasi-equilibrium, with the particles moving 
inside the cages formed by nearly frozen neighbors and 
the temperature measured from particles fluctuations is 
that of the heat bath. When t — t w > t w the system 
falls out of equilibrium, entering the aging regime, and 
the data in Fig. [I] depart from the FDT line. From the 
figure it is clear that the FDR is still a constant, but now 
< X(t,t w ) < 1. The constancy of the FDR in the out 



of equilibrium regime is one of the central predictions of 
MF approaches. Here we see that this is still valid for a 
finite dimensional system. 



N = 
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n C(t,t w ) 

FIG. 2. Same as Fig. hi, now with different volumes in order 
to show the absence of finite size effects. The linear fit to the 
data is very good and gives a value for the FDR X ~ 0.65. 
For the EA model (inset) the linear fit is far from good. 

In Fig. we compare the results for the FDR mea- 
sured on two large systems, whose sizes are N = 30 3 and 
N = 60 3 . Both curves correspond to t w = 10 4 and no 
finite size effects are evident. Fitting the data in the out 
of equilibrium regime to the straight line 

Tn(t,t w ) = X(t w )[q EA (t w ) - C(t,t w )] + [1 - q EA {t w )] , 

(11) 

allows us to compute the ^-dependent FDR X(t w ) and 
the Edwards- Anderson order parameter q EA (t w ). In the 
large times limit they should converge to the correspond- 
ing equilibrium values 

We report the linear fit in Fig. || in order to show how 
well the data can be fitted with the formula in Eq.(^Tj). 
As one can see from Fig. [j], the slope X(t w ) changes very 
little with t w and it takes the same value (within the er- 
ror) for the two largest waiting times. The results of our 
fits give X = 0.64(3) and q EA — 0.92(1). In comparison 
with other works the correlation range we are explor- 
ing may seem quite small. However it should be kept in 
mind that we are using non-connected correlations func- 
tions which tend in the large times limit to p 2 ~ 0.44. So 
we are actually spanning half of the allowed range. 

In the inset of Fig. || we present the same kind of data 
(t w = 10 5 , 10 6 ) for the Edwards-Anderson (EA) model, 
which is expected to have more than two time scales. It 
is clear that in the EA model the slope changes along 
the curve and a straight line is not able to fit the whole 
set of data. Moreover higher temperatures data [|ll| sug- 
gest that the slope still have to decrease (in modulus) for 
smaller correlations, making the linear fit even poorer. 

An important difference between the present model 
and MF approaches should be evident: activated pro- 
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cesses present in finite-dimensional systems should dom- 
inate the asymptotic dynamics and as t w — > oo equi- 
librium dynamics should be restored (a similar behavior 
can be observed in MF models by considering small sys- 
tems [^H ) . However the time a macroscopic system needs 
to reach such an equilibrium (thermalization time, i cq ) 
increases very rapidly with the system size and then the 
use of large sizes (here we have L = 30 and 60) prevents 
the system from reaching equilibrium in accessible time 
scales. In other words, we expect that the FDR explic- 
itly depends on t w and it tends to 1 in the limit t w — > oo. 
However in the range 1 <C t w <C t Qq (where we actually 
are) the FDR should relax into some very long plateau. 
It is very remarkable that finite times effects are very well 
described by MF theories. 



ratio X which is, with a very good approximation, in- 
dependent of temperature in the glass phase. This last 
observation does not agree with MF predictions. Being 
a model in 3D, the frustrated lattice gas is an interesting 
test ground for performing a systematic study of acti- 
vated processes, a main ingredient absent in MF models. 

This work was partly supported by Brazilian agen- 
cies CNPq and FAPEMIG. JJA acknowledges the Abdus 
Salam ICTP (Trieste) for support during his stay, where 
part of this work was done. 
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FIG. 3. Same as Fig. u\ for different temperatures in the 
glassy phase. 

In Fig. we show the usual integrated response ver- 
sus correlation plot for different temperatures, being the 
chemical potential always equal to fi = I. We estimated, 
both from the density measurements and from the FDR, 
the glassy transition to be located around T g ~ 0.2. So 
all the data refer to the glassy phase. The main result to 
be noted is the good parallelism between all the curves 
in the aging regime. They can be perfectly fitted with 
the same value for X and different values for q EA , that 
increases lowering the temperature. As T approaches T g 
from below, FDT is recovered. Also, as q EA (T) decreases, 
X remains constant within the numerical precision, what 
is in contrast with MF models where the FDR is usually 
proportional to the temperature (X oc T). The exact 
behavior near the transition, that is, whether X (and 
<7 EA ) are continuous or not |I6j, is difficult to establish 
numerically and will be addressed in a future work. 

In summary, we have found that the glassy phase of a 
realistic model of structural glass presents two well sepa- 
rated time scales, as found in MF models of spin glasses. 
The out of equilibrium long time dynamics can be charac- 
terized by a constant value of the fluctuation dissipation 
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